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We combine the concept of Bell measurements, in which two systems are projected into maxi- 
mally entangled states, with the concept of continuous measurements, which concerns the evolution 
of a continuously monitored quantum system. For time-continuous Bell measurements based on 
homodyne detection of light we derive the corresponding stochastic Schrodinger equations, as well 
as the unconditional feedback master equations. Our results cover in particular the two scenarios of 
time-continuous quantum teleportation and entanglement swapping. We apply our results to show 
that (i) two two-level systems can be deterministically entangled via homodyne detection of light, 
even including photon loss approaching the fundamental limit of 50%, and [ii) a quantum state of 
light can be continuously teleported to a mechanical oscillator. This time-continuous remote quan- 
tum state preparation works under the same conditions as are required for optomechanical ground 
state cooling. 



Introduction. — The quantum theory of time- 
continuous measurements provides profound insight 
into the intricate connection between decoherence. 
measurement and the transition from quantum to 
classical physics [1-12], see also [13-lG]. In conjunction 
with continuous feedback control it builds the basis 
for an increasing scope of quantum physics experi- 
ments, ranging from single atoms [17-19], cavity modes 
[20, 21[, via atomic ensembles [22-24[, up to massive 
mechanical oscillators [25-30[, demonstrating continuous 
measurement and control at the quantum limit. 

In such experiments the system to be measured and 
controlled (e.g., the ion, cavity, mechanical oscillator) 
continuously interacts with an environment serving as 
a 'meter system' (e.g., an optical field) which in turn is 
continuously measured. The corresponding measurement 
results are fed back to the primary system as a continu- 
ous control in order to achieve the desired quantum state 
or dynamics. The impressive experimental progress [17- 
30[ opens up the possibility to combine several of these 
setups, and to perform a continuous nonseparable mea- 
surement on the environments, i.e., the meter systems. 
Such a strategy was shown to be advantageous in the 
control of entanglement of parametrically interacting os- 
cillators [31, 32[. In the extreme limit such a measure- 
ment can correspond to a continuous projection of the 
environments on maximally entangled states, that is, to 
a continuous Bell measurement. 

In this article we introduce and study continuous Bell 
measurements realized via continuous homodyne detec- 
tion of light [5, G[. We derive the constitutive equations 
of motion — the conditional stochastic Schrodinger equa- 
tion and the unconditional feedback master equation — 
for the monitored systems. Two scenarios are studied 
(Fig. 1): In the first one a single system couples to a field 
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FIG. 1. Schematic setups for (a) time-continuous teleporta- 
tion and (b) entanglement swapping. The systems may take 
the form of (c) a harmonic oscillator (e. g., an optomechanical 
cavity), or (d) a two- level system (e.g., a single spin). 



A, which is subject to a Bell homodyne detection after 
interference with a second field B prepared in a general 
state of Gaussian (possibly squeezed) white noise. In 
the second setup the Bell measurement is performed on 
the fields emitted by two systems. Provided the system- 
field interactions create entanglement, the feedback mas- 
ter equations for the two setups cover in particular the 
cases of time-continuous quantum teleportation and time- 
continuous entanglement swapping. 

We illustrate the latter process for two two-level sys- 
tems, and show that they can be continuously and deter- 
ministically driven to an entangled state, ideally a Bell 
state. Entanglement is generated even for light losses up 
to 50% when the quantum capacity of the optical com- 
munication channel drops to zero [33[ . This deterministic 
scheme for distributing long distance entanglement can 
provide the basis for a dissipative quantum repeater ar- 
chitecture [34[. Furthermore, time-continuous quantum 
teleportation is exemplified for an optomechanical system 
where the quantum state of continuous-wave light is con- 
tinuously transferred to a moving mirror. We show that 
the experimental requirements are essentially the same 
as for optomechanical ground-state cooling [35, 36[, the 
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observation of back action noise [37], or ponderomotive 
squeezing ]38]. 

Continuous Teleportation. — We consider the setup 
shown in Fig. 1(a): A system S couples to a ID electro- 
magnetic field A via a linear interaction Hi^t — i[s 
s^a{t)], where s is a system operator, and the light field 
is described (in an interaction picture at a central fre- 
quency Wo) by a{t) = J doj a(a;) e"''-""""'*. Analogously 
a second ID field B is described by an operator b{t). Our 
first goal is to derive a stochastic master equation (SME) 
for the state of the system S, conditioned on the results 
of a time-continuous Bell measurement on the two fields 
A and B. In a Markov approximation we restrict our- 
selves to white noise fields, which means that both a{t) 
and b{t) are (5-correlated. This allows us to introduce the 
Ito increment dA{t) = a{t)dt (and analogously dA^ , dB, 
dSt), and to express the time evolution of the state |0) of 
the overall system (S+A+B) as a stochastic Schrodinger 
equation in Ito form [15, 39[, 

d|0) = (-ii/effdt + sdAt) 10), (1) 

where Hcs — Hsys — i-^s^s, with the (unspecified) system 
Hamiltonian TJgys- 

We assume that the initial state of the overall system 
is 10(0)) = |V'(0))s|vac)A|M)B, where |M) is an arbi- 
trary pure Gaussian state defined by the eigenvalue equa- 
tion [{N + M* + l)b{t) ~ {N + M)b^t)]\M)B = 0. The 
parameters N G R, M e C obey the relations > 
(A'' = M = true for vacuum) and |Mp = A^(A^ -I- 1). 
The present white noise model essentially assumes that 
the squeezing bandwidth is larger than all other system 
time scales. Making use of the fact that a{t)\(j){t)) — 
a{t)\(j){0)) = and the above eigenvalue equation, we 
can rewrite equation (1) in terms of the EPR operators 
X+ = {a+a'f+b+b'f)/V2a.ndP- = -i{a-a'' -b+b'f)/V2, 
which can be simultaneously measured in this setup. The 
resulting equation reads 

d|0) = [-ii/effdt + s i^idX+ + iiydP.)] 10), (2) 

with /i = (1 - M + M*)/{1 + N + M*) and v = 
i(l + 2N + M + M*)/{1 + N + M*). Writing Eq. (1) 
in this form enables us to project (2) onto the EPR state 
\I+I-)ab, defined by X+\I+I_)ab = I+\I+I-)ab and 
P_|/+/_)a_b ~ AB- This leads to the so-called 

linear stochastic Schrodinger equation [14[ 

d|i^c) = {-iifcffdi + s W+{t) + wl^{t)]dt} IV^c), (3) 

for the unnormalized system state I'i/'c), which is condi- 
tioned on the measurement results I± . Note that /_|_ and 
/_ are real- valued, Gaussian random processes, which are 
proportional to the measured homodyne photocurrents. 
As I± result from mixing the fields A and _B on a beam- 
splitter, they carry information about both fields, and 
will therefore be correlated, which is one of the crucial 
features of a Bell measurement. We can write [6, 15[ 

I+{t) = v/l72(s + st)^(^) + ^^(t), (4a) 
I.it)^iy/l/2{s~s^)^^t)+Ut), (4b) 



where (,±{t) = dW±{t)/dt is zero-mean, Gaussian, white 
noise with corresponding Wiener increments dW± [40[. 
The (co-) variances of dW± are given by 

wi dt:^{dW+f = [A^ + 1 + {M + M*)/2]dt, (5a) 

W2 dt-.^idW-)^ = [A + 1 - {M + M)/2]dt, (5b) 

W3 dt:^dW+dW- = -[i(M - M*)/2]di, (5c) 

as follows essentially from the initial mean values with 
respect to the optical fields ((A+)^)0(o), ((P_)^)0(o), 
etc. As expected, we in general find non-zero cross- 
correlations between /+ and /_, which depend on the 
squeezing properties of the input field B. Using Ito rules 
[40[ we can construct the corresponding stochastic master 
equation (in Ito form) for the system state conditioned 
on the Bell measurement result, 

dpc = Cpcdt+^{n[fis]p,dw++n[ii^s]p,dw^}, (6) 

where we defined Cp — — i[i/sys, p] + 2?[s]p, the Lindblad 
operator = (sps^ — ^ps^s — ^s^sp), and 'H[s]p = 

is-{s))p + p{s-{s))\ 

We now apply Hamiltonian feedback proportional to 
the homodyne photocurrents to the system, a scenario 
which covers in particular the case of continuous quan- 
tum teleportation of the state of field B to the system S. 
We follow [41[ in order to derive the corresponding uncon- 
ditional feedback master equation. Hamiltonian feedback 
is described by a term [pc]fb = + I-Kl-) Pc, 

where we define IC±p = —i[F±,p], and Hermitian oper- 
ators F±. After incorporating this feedback term into 
the SME (G) (see supplemental information), and taking 
the classical average over all possible measurement out- 
comes, p = E[pc\, we arrive at the unconditional feedback 
master equation 

p = -i [H^ys + (1/4) {(F+ +iF^)s + sHF+ - iF-)} , p] 

+ (l/2){l?[s - iF+]p + V[s - F^]p + W3V[F+ + F_]p 

+ {wi -ws- l)V[F+]p+{w2 ~W3~ l)V[F^]p}. (7) 

This is the main result of this section. The evolution 
of the system S thus effectively depends on the state of 
the field B (via Wi) which has never interacted with S, 
and which can in principle even change (adiabatically) in 
time. Eq. (7) can thus be viewed as a continuous "remote 
preparation" of quantum states. 

In order to illustrate this point we consider the case 
where the target system 5 is a bosonic mode. For a 
system to be amenable to continuous teleportation the 
system-field interaction must enable entanglement cre- 
ation. We thus set s — , with c a bosonic annihilation 
operator, and therefore obtain Hi^t oc ca{t) + c^a^{t), 
an interaction commonly known as two-mode-squeezing 
interaction. Additionally we choose F+ = i(c — c^) and 
F_ = {c+c'), which means that the photocurrents /+, /_ 
will be fed-back to the x and p quadrature, respectively. 
The resulting equation can brought into the form 

p = -i[ffsys, p] + {2N + l)2?[J]p, (8) 
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FIG. 2. (Color online) Performance of continuous entan- 
glement swapping: (a) Entanglement (logarithmic negativity 
[42]) of the stationary state versus z parameterizing the light- 
matter interaction oc \\{\J z{\ -|- z)(t^ — \/\ — 2cr~)a^(f) — h.c], 
with transmissivity r\ as indicated, (b) Entanglement versus z 
and r\ for optimized gains, (c) Entanglement versus transmis- 
sivity for optimized z. Nonzero entanglement can be achieved 
even for losses approaching 50%. 

where the jump operator J is determined by J cx — i(2iV+ 
1 — M — M* ) a; + ( 1 + M - M* )p (with an appropriate nor- 
malization). For Tfgys = equation ((S) has the steady- 
state solution yOss = |V')('0li where Jji/') = 0. Up to a triv- 
ial transformation this state is identical to the input state 
\M){M\. Note that for the vacuum case iV = M = we 
find J = c, which means that, devoid of other decoher- 
ence terms, the system will be driven to its ground state. 
Below, we will come back to this scenario, and discuss its 
implementation on the basis of an optomechanical system 
in more detail. First, however, we consider the second re- 
alization of quantum dynamics conditioned on a continu- 
ous Bell measurement, namely continuous entanglement 
swapping, as depicted in Fig. 1(b). 

Continuous Entanglement Swapping. — We now re- 
place the Gaussian input state in mode B with a field 
state emitted by a second system, which couples to the 
field B via i?int = i[s2b\t) — s\b{t)]. Using the same logic 
as before we can derive the linear stochastic Schrodinger 
equation for the bipartite state \'ip) (comprising 5*1 -I- ^2) 

d|V3) = [-ii^cffd^ + s+I+{t)dt + is_/_(i)dt] ivi), (9) 

where now Hcs = Hsyl + Hsyl — 5 J2i=i,2 ^l^i ^^'^ *± ^ 
Si ± S2- Accordingly, the homodyne currents read 

I+{t) = yi72(,s+ + (10a) 
I-{t) = 17172 (s_ - sl)^^t)+C-{t), (10b) 

and the corresponding conditional master equation is 

dpc = Cpcdt + ^/i/2 {n[s+]p, dw+ + n[is^]pc dw-} , 

(11) 

with Cp = ~i[H^yl + H^ylp] + V[si]p + V[s2]p. Here, 
the Wiener increments are uncorrelated and have unit 



variance, i.e., (dW+)2 = (dW-)^ = dt, dW+dW- = 0. 
Applying Hamiltonian feedback to either or both of the 
two systems in the same way as before gives rise to 

p^-i [H, p] - i(l/4) [{F+ + iF_)s, + h.c, p] 
-i(l/4) [(F+-iF_)s2+h.c.,p] 
+ (1/2) {V[s+ - iF+]p + V[s^ - . (12) 

This is the desired feedback master equation for continu- 
ous entanglement swapping. In a simple generalization of 
the continuous teleportation example considered above, 
one can check that the analogous strategy applied to two 
bosonic modes with Si = cj, will indeed drive the two 
systems to an Einstein-Podolsky-Rosen entangled sta- 
tionary state. In view of Fig. 1(b) the resulting topology 
comes close to a Michelson interferometer, for which a 
similar scheme was discussed in [43]. Note, however, that 
the central equations (6), (7) and (9), (12) are general, 
and also apply to non-Gaussian systems. As a rather 
surprising application we will show, that a pure entan- 
gled state of two two-level systems (TLS) can be created 
deterministically. 

Consider t wo TLS which couple to a ID field via oper- 
ators Si = y/ z{l + z)af — Vl — za^ {i = 1,2; z G [0, 1]). 
The fields are subject to a continuous Bell measure- 
ment as depicted in Fig. 1(b). The homodyne pho- 
tocurrents I±{t) are used in a Hamiltonian feedback 
scheme to generate rotations of the TLS about their 
X and y axes according to = G-crf + G+ct| and 
F_ = G+uf — G_cr|, with gain coefficients G± = 

y^z/{l + z) {l ± a/z(1 + z)/(l - z)) . For this choice of 
Si and F±, and assuming that the levels in each TLS 
are degenerate (such that i/iys = 0), the jump opera- 
tors in Eq. (12), become J+ = s+ — i-F+ oc ji — Aj2 and 
J_ = s_ — F_ (X j2 + Aji, where ji = + zct^ and 
j2 = + ^(^t ^ and A is a real coefficient. The common 
dark state of the jump operators J±\^) = ^1,2!^) = is 
the pure entangled |$) cx |00) — z|ll) which becomes a 
maximally entangled Bell state for z — 1 [34[. The par- 
ticular linear combination of ji^2 in J± is chosen such, 
that the state |$) is also an eigenstate of the effective 
Hamiltonian Hcff = \[{F+ + iF^)si + {F+ - iF_)s2 + 
h.c] - + j1j_] of Eq. (12), i.e., i7eff|$) = 0. 

Together, these properties guarantee that the stationary 
state of Eq. (12) is the pure entangled state 1$) [44[. Note 
that in this way entanglement is generated deterministi- 
cally, in contrast to conditional schemes based on photon 
counting [45-5 0[. Moreover, the present protocol does 
not require to couple nonclassical light, such as squeezed 
or single photon states, into cavities or to the TLS, in 
contrast to previous schemes [51-58[. 

The ideal limit of Bell state entanglement (z — > 1 [59[) 
is achieved only in the limit of infinite feedback gains 
G±, as is to be expected for the present treatment. More 
sophisticated descriptions of feedback, might relieve this 
restriction [15[. However, in the relevant case includ- 
ing losses, the optimal feedback gains stay finite even 
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in the present description. Assuming that all passive 
photon losses, such as finite transmission and detector 
efficiency, are combined in one transmissivity/efficiency 
parameter rj, we have to apply the generalized feedback 
master equation from the supplemental information in- 
stead of Eq. (12). For given rj and light matter inter- 
action, i.e., fixed z, we optimize the feedback gains G± 
in order to maximize the entanglement of the station- 
ary state pss- We keep the particular form of the feed- 
back Hamiltonians F± as it preserves the Bell diagonal 
structure of pss- Fig. 2 shows that entanglement can be 
achieved even for losses approaching 50%, which is where 
the quantum capacity of the lossy bosonic channel drops 
to zero [33]. 

The necessary strong coupling of TLS to a ID optical 
field can be achieved in a variety of physical systems, such 
as cavities [21, 60-65] or atomic ensembles [24, 66, 67]. 
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FIG. 3. (Color onHne) Mechanical squeezing ( against co- 
operativity C for varying mechanical bath occupation n = 
0, 1/10, 1/2, oo (represented by different colors) and sideband 
resolution n/ujm = 1 (10) [solid (dashed) lines]. The solid 
black line at (" = —6 dB shows the squeezing level of the input 
light, (corresponding to A'^ ~ 0.56). The dashed vertical line 
shows the value of Ccrit = l/[\/N[N + 1) - TV] « 2.7 above 
which mechanical squeezing is achievable for any n. 

Application to optomechanical systems. — In the re- 
mainder of this article we will show how continuous quan- 
tum teleportation can be implemented in an optomechan- 
ical system in the form of a strongly-driven Fabry-Perot 
cavity with one oscillating mirror [Fig. 1(c)]. Here the 
system Hamiltonian (in the laser frame at ajg) is ^fsys = 

^^0 + -ffom = (WniCj-nCm + Accjcc) -|- g{c^ + cl^){Cc + cj), 

where Wm is the mechanical frequency, Ac = Wc — wo is 
the detuning of the driving laser (at wq) with respect to 
the cavity (at Wc), and g is the optomechanical coupling 
strength. The decay rate of the cavity we denote by k, 
while the coupling of the mechanical oscillator to a ther- 
mal bath (with a mean occupation number n) we call 
7. Cm and Cc are bosonic annihilation operators of the 
mechanical and the optical mode respectively. 

In order to successfully implement continuous telepor- 
tation in this system we need to make several adaptations 



to the generic case presented above. To understand this 
we first need a clear picture of the system's dynamics: In 
the regime g <^ k <^ lo,^ and for A^. — — the optome- 
chanical interaction is i/om ~ g{cn\Cc + c^cj). Under the 
weak-coupling condition {g <^ k) the cavity follows the 
mechanical mode adiabatically, and we effectively obtain 
the required entangling interaction between the mirror 
and the outgoing field. Moreover, the mechanical oscil- 
lator resonantly scatters photons into the lower sideband 
at u!c- Spectrally, the photons which are correlated with 
the mechanical motion are therefore located at the side- 
band frequency iOc- Consequently, we have to modify 
the previous measurement setup in two ways: Firstly, we 
choose the center frequency of the squeezed input light 
to sit at the sideband frequency lOc- Secondly, we now 
use heterodyne detection to measure quadratures on the 
same sideband. These two modifications, together with 
the adiabatic elimination of the cavity (a perturbative 
expansion in g/n ]8]) and a rotating- wave approximation 
(an effective coarse-graining in time ]40]), allow us to 
write the SME for the mechanical system, in the rotat- 
ing frame at Wm, as 

dpc = 7-I?[ci„]/5cdt + 7+2?[4jpcdt (13) 
+ ^g^K/2 {H[-ipv+cl]PcdW++H[iyv+cl]PcdW^} , 

where r]± = [k/2 + i(Ac ± Wm)]^^- The first two terms 
describe passive cooling and heating effects via the op- 
tomechanical interaction with cooling and heating rates 
7_ = 7(n-|-l)-|-2g^Re(?7_) and 7+ = 7n-|-25^Re(77+), as 
was derived before in the quantum theory of optomechan- 
ical sideband cooling [68, 69]. The last two terms in (13) 
describe the continuous measurement in the sideband re- 
solved regime for arbitrary laser detuning Ac. This is 
an extension of the conditional master equation for op- 
tomechanical systems usually considered in the literature 
which concerns a resonant drive and the bad-cavity limit 
]10, 70] (see however [71]). 

For simplicity we assume here that we can apply feed- 
back directly to the mechanical oscillator. More re- 
alistic feedback, for example optical feedback through 
the cavity, will be considered elsewhere. Under this 
premise we can adopt the same choice of F± as before, 
and arrive at a feedback master equation similar to (8), 
P = j{n+l)V{c^]p+jnV[cl]p+{igyK){X,{e)V[Me)] + 
X2{eMJ2{e)]}p, where e = [1 + (4wm/K)2]-i. (For de- 
tails on how to derive and Ji refer to the supplemental 
information.) The protocol's performance is degraded 
by mechanical decoherence effects and counter-rotating 
terms of the optomechanical coupling, which are sup- 
pressed by e. For fixed input squeezing (determined by 
iV) it is possible to describe the dynamics in terms of 
the mechanical bath occupation n, the sideband reso- 
lution k/cJiji, and the so-called cooperativity parameter 
C ~ g'^/{n + 1)7^. To illustrate the performance of 
the protocol we plot the mechanical squeezing C, and 
compare it to the squeezing of the optical input state 
(Fig. 3). As is evident from the figure there exists a crit- 
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ical value Ccrit(A'') = l/[^yN{N + 1) - N] determined 
by the level of input squeezing, above which mechanical 
squeezing can be achieved for any thermal occupation n. 
We emphasize that this condition on the optomechanical 
cooperativity is essentially the same as for the recently 
observed ground-state cooling [35, 36], back-action noise 
[37], or ponderomotive squeezing [38]! 

Conclusion. — In this article we present a generaliza- 
tion of the standard continuous-variable Bell measure- 
ment based on homodyne detection to a continuous mea- 
surement setting. We show how this concept, together 
with continuous feedback, can be applied to extend ex- 
isting schemes for teleportation and entanglement swap- 
ping. We specifically apply these protocols to optome- 
chanical and two-level systems, and remarkably find, that 
in addition to EPR entanglement, we can also generate 
qubit entanglement in a deterministic fashion. We sug- 
gest that the formalism developed here can serve as a ba- 
sis for continuous measurement based quantum commu- 
nication and information processing with both discrete 
and continuous variables. 
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Appendix A: Bell measurement master equations 

We first treat the case of continuous teleportation. 
Note that due to their definition, the Ito increments 
commute with the unitary evolution operator at all 
times. It thus holds that dA{t)\(j){t)) = dA(t)|^(0)) = 
dA{t)\va.c)A = 0, and by the same reasoning [{N + M* + 
l)dB{t) - {N + M)dB'f {t)]\(f>{t)) = 0, for the initial state 
\<p{0)) = |'(/'(0))s|vac)yi|M)B- Inserting these terms into 
(1) with appropriate prefactors yields 

d|</>) = {-iiJeffdi + s(dA^ + adA + dB ~ adB'')} 

where a = {N + M)/{N + M* + 1). Rearranging this 
leads to Eq. (2). The probability distribution of the 
measurement results I± is given by T(J+(i), /_ (i)) = 
\{I+I-\cj){t + dt))\ ]6, 15], which, to first order in dt, is 
a Gaussian with first and second moments given by (4) 
and (5) respectively. 

To find the SME corresponding to we define Pc{t + 
dt) = \'il)c{t + dt)){i;cit + dt)\ and note that \ipcit + dt)) = 

\^c{t)) +d|'0c(O)i where {ipc) and are unnormalized. 
After normalizing pc{t + dt) = pc{t + dt)/tT[pc{t + dt)] 



we expand the resulting equation to second order in the 
noise increments dW± and apply the Ito rules (5) (see 
[40]). With the definition dpc{t) = Pc{t + dt) ~ p^{t) we 
find (()). 

We follow the procedure developed in ]41] to add the 
feedback term [pdfb = ■\/l/2 (/+/C+ + I-K. ) pc to the 
conditional master equation. Note that this term must 
be interpreted in the Stratonovich sense ]41]. To recon- 
cile it with equation (6) we thus have to convert (G) to 
Stratonovich form, add [pc]fb, and convert to result back 
to Ito form. This yields 

dpc = Cpcdt + \ [dX+/C+ + dP^K-f Pc 
+ \ [dA+/C+ + dP_/C_] {dW+nips] + dW-n{ws\}p^ 

+ ^ {dw+n{ps] + dw-n\vvs]} Pc 

+ ^[dA+/C++dP_/C_]pc, 

where the operator ordering KfH. was used in order to get 
a trace-preserving master equation ]41]. Using the fact 
that dX^dXj = dX^dW^ = dW.dWj together with (5) 
and taking the average with respect to the measurement 
outcomes p = E[pc] this equation can be brought into 
the form (7). 

The case for entanglement swapping can be treated 
analogously. For the full system we can write 

d|0) = [-ii/effd< + sidAt + s2dB^ 10) 

= [-iH.sdt + si {dA^ + dB) + S2 {dA + dB^)] |0), 

which, by projection onto the EPR basis, leads to (9). 
Note that the initial state here is assumed to be |0(O)) — 
\'>P{Q))siS2W^'^)aW3.c)b- The corresponding feedback 
equation is derived as for continuous teleportation, us- 
ing the multipHcation table (dW+)^ = (dW^)^ — dt, 
dT^+dT^_ = 0. 



Appendix B: Bell measurement for non-unit 
detector efficiency 

Passive losses due to inefficient detectors or imperfect 
transmission can be accounted for by introducing a com- 
bined transmissivity / efficiency < rj < 1 . The equations 
in the main text can be generalized in a straight-forward 
manner ]15]. In particular we find for the case of contin- 
uous teleportation the conditional master equation 

dpc = Cpcdt+^/^{H[ps]p^ dW++n[ws]pc dW_} , 
and corresponding photocurrents 

/+(t) = v/W2(.5 + st)^,(,) + e+(t), 

/_(0=iy^(s-st)^(,)-f^_(t). 
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Including feedback as [pcjfb = \/l/2r] (/+/C+ + /-/C-) Pc 
gives rise to the feedback master equation 

p = -i [H,y, + (1/4) + iF^)s + s\F+ - iF^)] , p] 
+ I {v[s - iF+]p + V[s - F_] + ^2?[F+ + F_] 

+ "^^1^^^ V[F+] + "^-^^"^ V[fAp. 
V V J 

Applying the same considerations to the case of entan- 
glement swapping leads to 

dp, = Cp.dt + ^v/^{n[s+]p, dW+ + H[is_]pc dW_} , 
and 

p=-i [H, p] - i(l/4) + iF_)si + h.c, p] 
-i(l/4) [{F+-iF_)s2 + h.c.,p] 
+ (1/2) - iF+]p + V[s_ - F_]p} 

+ {l~7j)/r^{V[F+]p + V[F_]p}, 

replacing Eq. (11) and (12) respectively. 



Appendix C: Diagonalization of non-Lindblad terms 

In general the feedback master equations (7) and (12) 
are not in Lindblad form as the prefactors of the opera- 
tors V can be negative. To cure this we can rewrite the 
non-unitary part of the evolution in terms of i? = {x,p)"^ 
as p = J2ij^'i3 [RiPR] - \pRjRi - \R]Rip), where A 
is a Hermitian matrix. By virtue of the eigenvalue de- 
composition of A we can write p = 'Ylii^iT^[Ji]p with 
Ji — Vi ■ R, where and (i = 1, 2) are the eigenvalues 
and eigenvectors of A respectively. 



Appendix D: Continuous optomechanical 
teleportation 

To derive Eq. (13) we start from the SME describ- 
ing heterodyne detection of the optomechanical cavity's 
output light. It can be obtained from equation (G) by 
replacing s — Cce'^'°*, where Ajo is the local oscillator 
detuning, and adding decoherence terms for the mechan- 
ical subsystem, which are due to its coupling to a ther- 
mal bath. After going into an interaction picture with 
Hq — Lu^cl^Cjn + Accjcc we can adiabatically eliminate 
the cavity mode, by expanding the equations in the small 
parameter g/n up to first order [8, 72]. Under this ap- 
proximation and after setting Ajo = ujm the SME takes 
the form 

dpc Ccpdt - cf [cin + cjjj, ypc - Pcv''] dt+ 
^5^2 {H[-ipe'"»*2/]pe dW+ + HK^^-^ylpc dW_ } 



where £niP = -i[uJmclnCni, p]+j{n+l)'D[cni]p+jnT>[cl^]p, 
and y = rj-Cm + ij+c^^. The steady-state mean amplitude 
of the intracavity field is (cc) = —ig{y). The heterodyne 
photocurrents can thus be obtained from (4) by the re- 
placement (s) — T' — i(7(?/)e"^'"*, thus 

/+ = -iygV2(2/e'"»* -h.c.)(i) +e+(i) 
/_ = V?^(ye'"»* -t-h.c.)(i) +e-(i) 

To apply the rotating wave approximation we go to the 
rotating frame of the mechanical oscillator (in the fol- 
lowing denoted by pc) [73] and take a time average over 
an interval St, which comprises many mechanical peri- 
ods, but is short on the system timescales in the rotating 
frame, i. e., 1/g ^ St ^ l/t^m- Under this assumption we 
can pull Pc out from under the integral and drop terms 
rotating at a frequency Wm. At the same time we treat St 
as infinitesimal on the system's timescale and therefore 
replace St ^ dt. This leaves us with the equation 

dpc = 7_I?[cm]pcdi -I- 7+2?[4jpcdt 
+ v/52«/2 {H[-ipv+cl]p,dWo+ +H[iyr,+cl]PcdWo-} . 

We also introduced the time-averaged Wiener increments 
dWo± — J dW±, which approximately fulfill (-I), as far 
as Pc is concerned. In principle this equation contains 
additional measurement terms corresponding to sideband 
modes at frequencies ±2win, which in a RWA are not 
correlated to the DC modes dWo± and were therefore 
neglected. By renaming dWo± dW± we arrive at (13). 
Note that by applying the RWA to the decoherence and 
the measurement terms consistently we assure that the 
resulting equation is a valid Belavkin equation ]74]. 

To do feedback we do the same coarse-graining proce- 
dure for the photocurrents /o± — J I±dt, and find 

/o+ = -iv/.g2«/2(?7+4 - h.c.){t) + ^+it), 
lo- = ^g^{v+cl + h.c.){t)+^_{t). 

Within these approximations the system is equivalent to 
the generic case in the main text. Applying the same 
feedback procedure we obtain 

Pm = {l{n + l)2?[c„i] 4- 7nl?[4i]} Pm 

+ {Ag^/n) {(1 + e)V[c^] + w^Vix,^ + p„] 

-I- {Wi -W3- l)X'[pm] + {W2 -W3- l)I?[Xin]}Pin, 

In view of the previous section we can diagonalize this 
equation to obtain 

{7(n+l)2?[c,J+7"2?[c|n]}p 
+ (4gV«){Ai(e)I?[Ji(e)]+A2(e)I?[J2(6)]}p, (Dl) 

where and are obtained from the eigenvalue decom- 
position of 

/.^^f u-2-i(l + e) -W3 + ^il + e)\ 
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where e = [1 + (4wm/K)^]"^- In the Hmit e ^- Eq. (Dl) 
reduces to (8), apart from the decoherence terms of the 
mechanical subsystem, which counteract the squeezing 
of the mechanical mode by driving it towards a thermal 



state. Operating the protocol in a regime of strong coop- 
erativity / n ^ ^(fi + 1) suppresses these perturbative 
effects. 
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